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ABSTRACT 

We present a study of a comparison of spin distributions of subhaloes found associated with 
a host halo. The subhaloes are found within two cosmological simulation families of MiUcy 
Way-like galaxies, namely the Aquarius and GHALO simulations. These two simulations use 
different gravity codes and cosmologies. We employ ten different substructure finders, which 
span a wide range of methodologies from simple overdensity in configuration space to full 6-d 
phase space analysis of particles. We subject the results to a common post-processing pipeline 
to analyse the results in a consistent manner, recovering the dimensionless spin parameter 
We find that spin distribution is an excellent indicator of how well the removal of background 
particles (unbinding) has been carried out. We also find that the spin distribution decreases 
for substructure the nearer they are to the host halo's, and that the value of the spin parameter 
rises with enclosed mass towards the edge of the substructure. Finally subhaloes are less 
rotationally supported than field haloes, with the peak of the spin distribution having a lower 
spin parameter. 

Key words: methods: A^-body simulations - galaxies: haloes - galaxies: evolution - cosmol- 
ogy: theory - dark matter 



1 INTRODUCTION 

Within the hierarchical galaxy formation model, dark matter haloes 
are thought to play the role of gravitational building blocks, 
within which baryonic diffuse matter collapses an d becomes de- 
tectable dWhite & Reeslll978l:lwhite & Fren]dll99lh . Gravitational 
processes that determine the abundance, the internal structure and 
kinematics, and the formation paths of these dark haloes within the 
cosmological framework, can be simulated in great detail using A'^- 



body methods. However, the condensation of gas associated with 
these haloes, eventually leading to stars and galaxies we see today, 
is still at the frontier of present research efforts. A first exploration 
of the (cosmological ) form ation of disc galaxies has been presented 
in lFall & Efstathioul ( Il980l) , where it was shown that galactic spin 
is linked to the surrounding larger scale structure (e.g. the parent 
halo). In particular, the general theory put forward by Fall & Ef- 
stathiou reproduces galactic discs with roughly the right sizes, if 
specific angular momentum is cons erved, as bary ons contract to 
form a disc (previously suggested bv lMestell l ll963h ) and if baryons 
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and dark matter initially share the same distribution of specific an- 
gular momentum. 

While the theory has subsequently been refined, it always 
included (and still includes) such a coupling between the par- 
ent halo's angular momentum and the resulting galactic disc 
(cf.[Dalcanton et al.Vl997):'Mo et al.'('l998^;'N avarro & Steinmet3 
( I2OOO ): Abadiet ah (2003); Bett et al. (2010,) ). The origin of 
the halo's spin can now be understood in terms of tidal torque 
theory in which protohaloes g ain angu la r mom e ntum f r om th e 
surrounding shear field ( e.g., IPeeblesI l ll969l) : I White! ( Il984h : 
iBames & Efstathioul ( Il987h ) as well as by the build-up of angu- 
lar momentum through the c umulative transfer of angular momen- 
tum from subhalo accretion l lVitvitska et al J2002h . Whichever way 
the halo gains its spin, it is a crucial ingredient for g alaxy for- 
mation and all semi-analytical modelling of it ( Kauffm ann et al.l 



1993','l997';'Frenk et al.'l997'; Cole et al. 2000; Benson et al. 2001 
Croton et al...2006; .De Lucia & Blaizot, ,2007.; .Bower et al.ii200€ 



Bertone et al.ll2007l ; lFont et alJl2008l ; lBensorJl2012l) 



A number of studies have been performed o n the spin of 
haloes , i n particular studies by IPeeblej ( Il969h; Bullock et al.l 
(l200lh; 'Hetznecker& Burkerj OOOeh; lBettetal.1 (l2007h; 
iMaccio et al. ( 2007); Gottlober & YepeT(l2007h; iKnebe & PoweJ 



(l2008h; lAntotiuccio-Delogu et al.l j2010l)~ IWang et al 
iTrowland et all j2012l) ; iLacerna & Padillj ilQm 



(2011); 



Brvan et aP 

( l2012h but so far little has been done on subhaloes. These studies 
look at the spin of individual dark matter haloes found in cosmo- 
logical simulations and generally do not focus on the substructure, 
or differences between substructure definition due to lack of 
resolution. Here we present a comparison of spin parameters 
across a number of detected subhaloes found by a variety of 
substructure finders. The finders use many different techniques to 
detect substructure within a larger host halo. This is a follow-up 
to a more general pa per comparing the r ecovery of structure 
by different finders in lOnions et all j2012h and its predecessor 
lRnebeetal.l ( l201ll) . 

The techniques studied here for finding substructures include 
real-space, phase-space, velocity-space finders, as well as finders 
employing a Voronoi tessellation, tracking haloes across time us- 
ing snapshots, friends-of-friends techniques, and refined meshes as 
the starting point for locating substructure. With such a variety of 
mechanisms and algorithms, there is little chance of any systematic 
source of errors in the collection of substructure distorting the re- 
sult. Subhaloes are particularly subject to distortion and evolution, 
more so than haloes because, by definition, they reside within a host 
halo with which they tidally interact. This can affect their structure 
and other parameters, and in this case we are particularly interested 
in the spin properties. We quantify the spin with the parameter A, 
a dimensionless quantity that characterises the spin properties of a 
halo and is explained in more detail in Section[2] 

The rest of the paper is structured as follows. We first describe 
the methods used to quantify the spin of the halo in Section|2] The 
data we used is described in Section[3] Next we look at the overall 
properties of the spin in Section|4T| Then we look at the correlation 
between the host halo and the subhaloes spin in subsection 14.21 
Finally we look at how the spin is built up within the subhalo as a 
function of mass in subsection l4.3l We conclude in Section|5] 



2 METHOD 

2.1 Spin parameter 

The dimensionless spin parameter gives an indication of how much 
a gravitationally bound collection of particles is supported in equi- 
librium via net rotation compared to its internal velocity dispersion. 
The spin parameter varies between 0, for a structure negligibly sup- 
ported by rotation, to values of order 1 where it is completely ro- 
tationally supported, and in pra ctice maximum values are usually 
A ^ 0.4 dFrenk & Whitell2012h . Values larger than 1 are unstable 
structures not in equilibrium. 

The re are two vari ants of the spin parameter that are in com- 
mon use. |Peeblesl l [l969l) proposed to parametrise the spin using the 
expression given in Equation[T] 



A : 



J^/\E\ 



(1) 



where J is total angular momentum, E the energy and M the mass 
of the structure. In isolated haloes, all of these quantities are con- 
served, which g ives th e definition a time independence. 



iBett et"aLril2007l) measured the Peeble's spin parameter and 
fitted an expression to th e distribution for hal oes extracted from the 
Millennium simulation dSpringel et al.ll2005h ; that is characterised 
by Equation [2] 



P(log ^) = ^ ( 1^ ) exp 



where A is 



An. 



yl = 31nl0 



r(a) 



(2) 



(3) 



The variables Ao and a are free parameters, and r(Q) is the gamma 
function. The best fit they found for field haloes was with Ao = 
0.04 326 and a = 2.509 

iBullock et al. I( l200ll) proposed a different definition of the spin 
parameter. A', expressed in Equation [4] As it is not dependent on 
measuring the energy it is somewhat faster to calculate when deal- 
ing with large numbers of haloes. 



A' 



J 



V2MRV 



(4) 



Here J is the angular momentum within the enclosing sphere of 
virial radius R and virial mass M, and V is the circular velocity 
at the virial radius {V^ — GAI/R). The Bullock spin parameter 
is more robust to the position of the outer radius of the structure. 
Bullock proposes a fitting function to the dist ribution as described 
in Equ ation \5\ which was based on one from lBames & Efstathiotj 
( Il987h . 



P(A') = 



AV27rcr 



exp 



ln^(A7A;.) 
2cr2 



(5) 



This has free parameters Aq and a and lBullock et al.l ( l200lh found 
a best fit for field haloes at values of Aq = 0.035 and a — 0.5. 

The Peebles calculation is perhaps more well defined for a 
given set of particles, as it is calculated directly from the parti- 
cles properties, whereas the Bullock parameter is easier to calcu- 
late from gross halo statistics, and is not dependant on the density 
profile. Fo r more comparisons of the tw o parameters the reader is 
referred to lHetznecker & Burker3 j2006l) 
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2.2 The SubHalo Finders 

In this section we briefly list the halo finders that took part in the 
comparison project. More detai ls about the specific algorithms are 
available in lOnions etal] ( l2012h and the articles referenced therein. 

A DAPTAHOP (Tweed) is a configuratio n space over density 



finder jAubert 



lOP (Tweed ) 
etalj2004l :n 



ailbopgh . 

• AHF (Knollmann & Knebe) is a con figuration space 
spherical overd ensity adaptive mesh finder jGill et all |2004| : 
iKnoUmann & K nebe 200^). 

• GRASSHOPPER (GRadient ASSisted HO P) (Stadel) is a re- 
working of the SKID group finderl lStadeteoOll) and appears within 
our wider comparison for the first time here, and so is described 
in more detail. It take s an approach like the HOP algorithm 
teisenstein & Hull Il998l) and reproduces the initial grouping of 
SKID in two computational steps. First densities are calculated for 
all particles as before using the Monaghan M3 SPH kernel over 80 
nearest neighbours. Second, for each particle, the gradient of the 
density is calc ulated in a way th at cancels the so-called EO error 
in the gradient jRead et alj|2010l) . by using the gradient of the M3 
kernel. Then, given that the neighbours lie within a ball of radius 
2/i, we create a link from this particle to the closest neighbour to 
the point a distance h in the direction of the gradient. 

After links have been created, each particle follows the chain of 
links until is reaches a cycle, marking oscillation about the den- 
sity peak of the group. Finally since noise below a gravitational 
softening length causes a lot of artificial density peaks we search 
for particles of a cycle which are within a distance r of any other 
particles in a cycle. The parameter r is typically set to 4 times the 
gravitational softening length, as was the typical case for SKID. Un- 
binding is also performed in a nearly equivalent way to SKID, but 
now scales as 0{n log n) as opposed to 0{n^) as was the case with 
the original SKID. 

The group finding with GRASSHOPPER is now fast enough to al- 
low it to be performed during a simulation but gives nearly identical 
results to the previous SKID algorithm. 

• Hierarchical Bound-Tracing (HBT) (Han) is a tracking algo- 
rithm working in the time domain jHan et al ]|201lh . 

• HOT-l-FiEstAS (H0T3d & H0T6d) (Ascasibar) is a general- 
purpose clus tering analysis tool , working either in config uration or 
phase space l lAscasibar & Binn ev 2005; Ascasibarll2010h . 

• MENDIETA (Sgro, Ruiz & Merchan) is a Friends-of-Frien ds 
based finder that works in configuration space jSgro et al.l2010l) . 

• ROCKSTAR (B ehroozi) is a phase-space halo finder 
teehroozi et~al]|201 ll) . 

• STF (Elahi) is a velocity space/phase-space finder jElahi et al.l 

limh . 

• SUBFIND (Spr ingel) is a configuration space finder 
dSpringel et aLll200lh . 

• VOBOZ (Neyrin ck) is a Voronoi tessellation based finder 
dNevrinck et alJ2005h . 



Tweed et 



3 THE DATA 

3.1 Simulation Data 

The firs t data set used in th is paper forms part of the Aquarius 
project l lSpringel et al.ll200^ . It consists of multiple dark matter 
only re-simulations of a Milky Way-like halo at a variet y of reso- 
lution s performed using GADGETS (based on G ADGET2. ISpringell 
l2005l) . We have used in the main the Aquarius- A to E halo dataset at 



Table 1. Summary of the key numbers in the Aquarius and GHALO simula- 
tions used in this study. A^high is the number of particles with the highest 
resolution (lowest individual mass). N250 is the number of high resolution 
particles found within a sphere of radius 250 kpc/h from the fiducial centre 
at each resolution {i.e. those of primary interest for this study). 



Simulation 


^high 


N250 


Aq-A-5 


2,316,893 


712,232 


Aq-A-4 


18,535,97 


5,715,467 


Aq-A-3 


148,285,000 


45,150,166 


Aq-A-2 


531,570,000 


162,527,280 


Aq-A-1 


4,252,607,000 


1,306,256,871 


Aq-B-4 


18,949,101 


4,771,239 


Aq-C-4 


26,679,146 


6,423,136 


Aq-D-4 


20,455,156 


8,327,811 


Aq-E-4 


17,159,996 


5,819,864 


GH-4 


11,254,149 


1,723,372 


GH-3 


141,232,695 


47,005,813 



2 = for this project. This provides 5 levels of resol ution, varying 
in complexity for which further details are available in lOnions et al.l 
(12012) . 

The underlying cosmology for the Aquarius si mulations is the 
same as that used for the Millennium simulation dSoringel et al.l 
I2OO5I) i.e. Q.M = 0.25, = 0.75, as = 0.9, n, = 1, h = 0.73 
These parameters are close to the latest WMAP data jjarosik et al.l 
[20111) (D.M = 0.2669, = 0.734, erg = 0.801, = 
0.963, h = 0.71) although erg is a little high. All the simulations 
were started at an initial redshift of 127. Precise det ails on the set- 
up and performance of these models can be found in lSpringel et al.l 
( l2008h . 

The second d ata set was from the GHALO simulation data 
jStadel et alj l2009l) . GHALO uses a slightly different cosmol- 
ogy to Aquarius, Qm = 0.237, Qa ~ 0.763, erg — 
0.742, Us = 0.951, h — 0.735 which again are reasonably close 
to WM AP latest results. It also uses a different gravity solver, PKD- 
GRAV2( [stadel et al.l2002h . to run the simulation therefore allowing 
comparison which is independent of gravity solver and to some ex- 
tent the exact cosmology. 

The details of both simulations are summarised in Table [T] 



3.2 Post-processing pipeline 

The participants were asked to run their subhalo finders on the sup- 
plied data and to return a catalogue listing the substructures they 
found. Specifically they were asked to return a list of uniquely 
identified substructures together with a list of all particles associ- 
ated with each subhalo. The broad statistics of the haloes found are 
summarised in Table|2] 

To enable a direct comparison, all the data retur ned was sub- 
ject to a common post-processing pipeline detailed in lOnions et al.l 
( l2012h . For this project we added a common unbinding procedure 
based on the algorithm from the AHF finder which is based on 
spherical unbinding from the centre. We requested data to be re- 
turned both with and without unbinding to allow a comparison of 
that procedure to feature in this study. Unbinding is the process 
where the collection of gathered particles is examined to discard 
those which are not gravitationally bound to the structure. This 
common unbinding allowed us to remove some of the sources of 
scatter introduced by the finders using slightly different algorithms 
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Table 2. The number of subhaloes containing 300 or more particles and centres within a sphere of radius 250kpc/h from the fiducial centre found by each 
finder after standardised post-processing (see Section lX2t . 



Name 


ADAPTAHOP 


AHF 


GRASSHOPPER 


HBT 


HOT3D 


H0T6D 


MENDIETA 


ROCKSTAR 


STF 


SUBFIND 


VOBOZ 


Aq-A-5 


24 


23 


23 


23 


18 


23 


17 


25 


22 


23 


21 


Aq-A-4 


222 


189 


170 


169 


174 


176 


123 


182 


155 


154 


163 


Aq-A-3 




1259 


1202 


1217 






787 


1252 


1124 


1117 


1141 


Aq-A-2 




4230 




4036 








4161 




3661 




Aq-A-1 




30694 












25009 




26155 




Aq-B-4 




197 




191 








202 




188 




Aq-C-4 




152 




146 








158 




137 




Aq-D-4 




217 




216 








230 




196 




Aq-E-4 




218 




219 








221 




205 




GH-4 




58 


58 










60 


54 


54 




GH-3 




1172 


1148 










1148 


1033 


1090 





— AdaptaHOP 


H6D 


STF 


— AHF 


HBT 


— subfind 


— Grasshopper 


— mendieta 


VOBOZ 


— H3D 


Rockstar 
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Figure 1. A comparison of the Peebles and Bullock spin parameters against 
t'max based on all finders using a common unbinding procedure from sub- 
haloes with more than 300 particles. The mean value of A/ A' is shown to- 
gether with one standard deviation error bars. It shows there is a correlation 
between the two but not a one-to-one correspondence, with some scatter 
present. The scatter at low t)max where haloes have very few particles is 
particularly pronounced. 



for removing unbound particles and to find what difference this 
made to the results. 

Both the halo finder catalogues (alongside the particle ID lists) 
and our post-processing software are available from the authors on 
request. 



4 RESULTS 

The results used were restricted to subhaloes with more than 300 
particles, as these produce a relatively stable value for s pin. Values 
below this limit tend not to converge across resolutions teett et al] 
l2007h . 





















/ 





X/kpc/h 





x/kpc/h 



Figure 2. An example of the influence of unbinding. Left panel: particles 
in the object prior to unbinding. Right panel: particles in the object after 
unbinding the been performed. The vectors indicate the direction and ve- 
locity relative to the bulk velocity of the individual particles making up this 
example subhalo. The contribution from the background particles has only 
a minor influence on the mass and Vmax of the subhalo, but a large effect 
on the spin parameter. 



4.1 Spin parameter 

In general there is a proportional relationship between the Pee- 
bles and Bullock spin parameters recovered by all the find- 
ers for the same subhaloes, although there is some scatter as 
shown in Figure [T] We do not dwell on the differences between 
the two definitions as that has already been studied elsewhere 
l lHetznecker & Burkertll200^ . As both definitions of spin exist in 
the literature we consider both metrics when comparing how the 
spin is recovered across finders, placing particular emphasis on 
their application to subhaloes. 

The majority of field haloes are found to c luster around a v alue 
of Ao — 0.044 for the Peebles spin para meter (Be lt et al. 20071) and 
Aq = 0.035 for the Bullock parameter l lBullock et al.a2001ll with a 
spread of values matched by a free parameter to give the width of 
the distribution. 



4.1.1 Spin for subhaloes with no unbinding performed 

If unbinding has not been correctly implemented the high speed 
background particles can distort the spin parameter enormously. 

To emphasise the type of structures that are found, an example 
of a subhalo without (left panel) and with (right panel) unbinding is 
shown in Figure[2] This is displayed as a vector plot of all the com- 
ponent particles position and velocities that make up the subhalo 
with the velocity vectors scaled in the same way in both panels. 
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The bulk velocity of the subhalo has been removed and all posi- 
tions and velocities are relative to the rest frame of the subhalo. 
Evident in the left panel of Figure |2] without unbinding are stray 
particles that are part of the background halo. Despite their small 
number these particles have both a large lever arm and large veloc- 
ity relative to the halo, and significantly alter the derived value of 
the spin parameter due to their large angular momenta. 

Comparing the two forms of the spin parameter in Figure [3] 
and Figure |4] we show how the spin parameter is quite chaotic, 
not matching a smooth Gaussian like profile as might be expected, 
and is clearly a long way removed from the idealised curve others 
have found for the distribution of spin. A significant number of the 
haloes have spin parameter values above 1, which is unphysical as 
these objects would be ripped apart by this level of rotation and so 
clearly cannot be equilibrium systems. This result is perhaps not 
surprising given the contribution from unbound background parti- 
cles moving with velocities far from the mean of the object being 
considered but clearly shows how poor unbinding methods are rel- 
atively easy to detect by looking at the spin parameter distribution. 
The Peebles spin parameter is more affected by the lack of unbind- 
ing than the equivalent Bullock parameter as it takes into account 
the kinetic energy of all the particles. Some more objective num- 
bers for this and subsequent comparisons are given in Table[3] 

The best fit values shown by the bold dashed lines are vastly 
different from the fiducial values given in Section [J] It is however 
significant that the finders HOT6d, ROCKSTAR and STF (shown by 
dotted lines) which all have a phase space based component in 
their particle collection algorithm already show a much better fit 
to the fiducial value than the non phase-space finders. It should be 
noted that when GRASSHOPPER is run without unbinding, it finds 
a large number of subhaloes which would normally be discarded 
by the unbinding procedure that is integral to the final part of the 
GRASSHOPPER algorithm. 



4.1.2 Spin for subhaloes with finders own unbinding performed 

Including each finder's own unbinding procedure improves the spin 
parameter measure considerably, as shown in Figure [5] and Fig- 
ure |6] Note that as ADAPTAHOP doesn't do any unbinding in its 
post-processing steps it is a clear outlier on this plot. The MENDI- 
ETA finder shows a double peak, which is indicative of some of the 
unbinding failing, an issue that the authors of the finder are cur- 
rently working on. 

When fitting the best fit curves to this data obtained for the 
spin parameter of subhaloes, the peak of the Bullock fitting curve 
given in Equation |4] is less than the field halo value by about 20 
percent, offsetting the mean towards smaller values of the spin pa- 
rameter. For the Peebles spin parameter the best fit is again offset by 
about 36 percent from the field halo value, again towards a smaller 
value of the spin parameter. 



4.1.3 Spin for subhaloes with a common unbinding performed 

Once a common unbinding is done, the curves move significantly 
closer to the idealised curve, although there is still some separa- 
tion. The plots of Figure |7] and Figure[8]compare the spin parame- 
ter distribution of the different finders using a common unbinding 
process. It shows th e m atch between the best fit curve quoted in 
iBullock et aD bOOlh and lSett et al] bOOTh and the haloes found by 
the finders taking part in the comparison. The values are now offset 
by 10 percent for the Bullock fit, and 30 percent for the Peebles 

© 2012 RAS, MNRAS OOOQ.mflOl 



Table 3. Summary of the best fit parameters for the graphs shown. Shown 
are the values for Aq and the other free pai'ameter (a or cr) used in the 
best fit, and their difference from the published field halo fit value. The 
subscripts F, N, O and C are for field haloes, no unbinding, own unbinding 
and common unbinding respectively. The A values are the difference from 
the field halo values, and the change is the percentage difference. All results 
are for level 4 data except the last which is level 1 



Plot 


Ao 


AAo 


change 


a /a 


Afr/a 


change 


Bullock J? 


0.035 






0.5 






PeebleSi? 


0.044 






2.509 






Bullock jv 


1.646 


1.611 


+4600% 


1.36 


0.86 


172% 


Peebles jv 


12.6 


12.573 


+29000% 


41 


39.2 


1560% 


Bullocko 


0.028 


-0.007 


-20% 


0.727 


0.227 


45.5% 


Peebles o 


0.028 


-0.016 


-36% 


3.643 


1.134 


45.2% 


Bullockc 


0.031 


-0.004 


-10.4% 


0.75 


0.25 


50.0% 


Peebles (7 


0.03 


-0.013 


-30% 


3.96 


1.448 


57.7% 


Bullock-Llo 


0.022 


-0.013 


-38% 


0.693 


0.193 


38.6% 



fit. This results in the closest fit to the data, although the subhalo 
spin again extends to slightly lower values for both parameters, and 
follows the best fit line at larger values. These results also have a 
similar trend for the Aquarius B-E haloes and the GHALO data sets. 
These inclusions show that the results are not influenced greatly by 
the simulation, simulation engine or small changes in the cosmol- 
ogy used. 

4.1.4 Spin at higher resolutions 

Going to higher resolutions afforded by the level I data as shown 
in Figure |9] the trend to a lower spin distribution peak continues, 
although only three of the finders were able to manage such a com- 
putationally intensive task. 

There is a more pronounced tendency to depart from the field 
halo fit line at low spin part of the distribution, with the peak and 
bulk of the distribution moving towards lower spin parameter val- 
ues. The finders also show more scatter with each of them identi- 
fying the peak of the distribution in slightly different places. The 
agreement particularly at the low end of the spin distribution is 
good but with slightly lesser agreement at the high end. 

Although AHF appears to find slightly more higher spin 
haloes, this is a result of the spherical unbinding algorithm it uses, 
which tends to also increase the spin distribution of the other find- 
ers slightly when used as the common unbinding procedure. 

The dashed line representing the level 4 data is included to 
allow a direct comparison between the level 4 and level 1 average 
fits. It shows the continued movement of the distribution towards 
lower spin values with higher resolution and an increase in data. 

4.1.5 Spin distribution summary 

The best fit curve figures for all these plots are summarised in Ta- 
ble [3] Even after cleaning the catalogues significantly by utilising 
a common unbinding procedure for all finders there remains a def- 
inite trend for substructure spin to be less than that found for field 
haloes. We investigate the reason for this in the next sections. 

4.2 Host halo radial comparison 

Next we consider whether the location of a subhalo within a host 
halo has any effect on the recovered spin parameter. First we 
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Figure 5. The same plot as Figure ^ but with the finders own unbinding 
processing applied to the data. This groups the spin parameters somewhat 
more tightly, and shows that spin is a good indicator of how well the un- 
binding procedure is removing spurious background particles. The ADAP- 
TAHOP finder doesn't perform an unbinding step, and this plot also shows 
up a flaw in MENDIETA's un binding procedure. Th e dashed line is the Bul- 
lock field halo fit curve from lBullock et al.l j200lh . The Bullock data fit is 
the best fit to the average using the Bullock fitting formula. 
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Figure 6. The same plot as Fi gure[5]except that this time the dashed line is 
the Peebles field halo fit from lBett etaP j2007h . The Peebles best data fit 
is the best fit to the average of the Bett formula. 



demonstrate in Figure [TO] that any effect is not an artefact of the 
finding process. Substructures closer into the centre of the host halo 
are more difficult to detect particularly by some finders, and there- 
fore subject to a loss of constituent particles that could be attached 
to the subhalo as shown in Muldrew et al. (201 10. To test this sup- 
position we took a subhalo found in the outskirts of the Aquarius-A 
main halo, and repositioning it at points closer to the location of the 
centre of the halo. Then two of the finders (AHF and ROCKSTAR) 
were rerun on the new data and the spin value calculated anew. The 
results shown in Figure[TO]indicate that there is little change in the 



value of the spin parameter with radius despite some variation in 
the recovered number of particles. 

Next we look at whether the mean value of the measured spin 
parameters changes with respect to the distance from the centre 
of the host halo. Figure [TT] displays this radial dependence for the 
indicated finders after a common unbinding step has been applied. 
The background points indicate the scatter in the spin parameter for 
any individual halo, as seen in the previous section. This shows a 
small trend for a lower mean spin as the subhaloes get closer to the 
centre of the host halo. This confirms the result that were found in 
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iReed et all ( [iooj) but is shown here at higher resolution and across 
more finders than the earher paper. 

Equivalent results are found when we compare 6 different sim- 
ulations generated by two different N-body codes and aggregate the 
average of the differen t finders across m ultiple haloes in Figure[T2l 
This effect (as noted in lReed et al .120041) is difficult to detect obser- 
vationally, as most substructure will form galaxies before falling in 
so will have its spin detect able from observations of galactic ro- 
tation curve already fixed ( lKauffmannetalJll993l) . The possible 
exception to this are galaxies forming at high redshift where the in- 
falling su bstructure has not yet formed stars, such as gas-rich dark 
galaxies jCantalupo et al.ll2012l) . made entirely of dark matter and 
gas, which may form structure after falling into a parent halo. 



4.3 Build up of the spin parameter within a subhalo 

This leads to the question of what causes the drop in the measured 
spin parameter with proximity to the centre of the host halo. Fig- 
ure [T3] shows the average change in the measured spin parameter 
as the detected subhalo is analysed from the centre outwards to its 
radius. This procedure is computed after the common processing 
and unbinding steps have been done. The subhaloes analysed in 
this way are then further binned into radial bins determined from 
the centre of the host halo. The outermost subhaloes, which are 
the least disrupted, show an initial decrease in measured spin pa- 
rameter as particles are removed from their outer edges. Subhaloes 
extracted from nearer the centre of the host halo do not show this 
initial decrease but instead have a monotonically rising spin param- 
eter as material is removed. 

This trend suggest that subhaloes are preferentially stripped of 
high angular momentum particles which are likely to be the most 



© 2012 RAS, MNRAS OOOQ.mflOl 



8 



Onions et al. 



0.07 
0.06 
0.05 
0.04 
0.03 
0.02 



AHF 

Grasshopper . 

H3D ■ ' 

H6D 

.HBT 



Rockstar 
STF 

sublind 
VOBOZ 




200 300 400 

Radial distance/kpc/h 



Figure 11. Comparison of mean spin parameter against radius from tlie 
centre of tlie host lialo. Common unbinding was applied in the pipeline in 
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above 300 particles are found there. The background points indicate the 
measured spin parameter for individual subhaloes. 
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Figure 13. The radial profile of the spin parameter across the the subhalo. 
This shows the change in the measured spin parameter as spin is analysed 
from the centre to the radius of the subhalo. Here Rmax is the subhaloes 
maximum radius. Each line represents a different host halo radial bin. Sub- 
haloes near the centre of the host halo show monotonically rising spin pa- 
rameter values spin, whereas further out the spin parameter initially drops 
before rising. 
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Figure 12. Comparison of mean spin parameter against radius from the 
centre of the host halo for several different haloes. The finders own un- 
binding procedure was used in the pipeline in this case. Each line is the 
average of the spin parameter binned into 10 bins across all finders partak- 
ing (AHF, GRASSHOPPER, ROCKSTAR, SUBFINDand STF). The haloes used 
were the Aquarius-A to E and GHALO all at level 4 of the resolution. The 
dashed/dotted lines indicate 20 and 80 maximum percentiles across all data. 



weakly bound particles, leading to a decrease in the spin parameter 
as they enter the host halo. The outermost particles are usually those 
least bound so are the most likely to be removed on infall. 

We can also examine how the spin parameter is built up as 
mass is added to a subhalo. In Figure[T4]we look at how the spin pa- 
rameter changes at various mass cuts of the subhalo, A/(< Mtot)- 
This shows how the spin is built up across the structure of the sub- 
halo. For each halo we calculate the spin parameter at 0.25, 0.5, 
0.75, 0.95 of the subhalo's total mass for all the contributing halo 
finders. We plot the mean and the standard deviation at each mass 
cut. 

As expected from Figure [O] all finders agree that the calcu- 
lated spin increases as the fraction of the subhalo mass that is used 




Figure 14. Comparison of the normalised mean Peebles spin at different 
mass shells of all subhaloes. The cuts were taken at 0.25, 0.5, 0.75, 0.95 
and the complete mass of the subhalo. A common unbinding procedure was 
run on the results. There is a clear decrease in spin with increasing contained 
mass, and about a 3-fold drop is evident. The top plot shows the value of 
the spin parameter, and the bottom plot the spin parameter normalised to the 
value of A at the subhaloes Rmax - Error bars are one standard deviation. 



to calculate the spin parameter is reduced. Note that haloes have 
steeply rising density profiles and so the inner 50% of the mass 
is contained within a much smaller fraction of the radius and that 
this result is averaged over all the recovered haloes and not split in 
radial bins. 
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5 SUMMARY & CONCLUSIONS 

There is a good level of agreement amongst the finders on the re- 
covery of the distribution of the spin of subhaloes, although differ- 
ences are still evident, causing scatter in some of the comparisons. 
Undoubtedly some of the scatter is due to different types of subhalo 
that are being recovered by the finders, some finders focusing on 
stream like structures and some on simple overdensities. There is 
still some room for improvement of the finders as the common un- 
binding test shows. Some of the possible improvements and sources 
of error will be outlined in Knebe et al. (in prep). 

The distribution of spin provides a very good indicator of the 
finders unbinding ability and seems broadly unaffected by the cos- 
mology and simulation engine in use. As such, the spin distribution 
serves as a mechanism to detect if substructure finders are perform- 
ing the unbinding correctly. The unbinding errors can be masked 
in other comparisons such as Wmax and mass plots but show up 
in an obvious way when the spin distribution is examined. Phase- 
space finders are less sensitive to poor unbinding as they have 
some implicit unbinding in the ir selection criteria when loo king 
at velocity com ponents. Indeed iHetznecker & Burker } i l2006h and 
iD'Onghia &">J avarro (20o3) both show there is a good correlation 
between the virialisation of haloes and the spin parameter, thus in- 
dicating its use for the determination of how relaxed the halo is, 
which is not unrelated to the unbinding process. 

The mean spin parameter of subhaloes decreases as they ap- 
proach the host halo's centre. This is a real effect and not an artefact 
of any difficulty in recovering structure as the subhalo approaches 
the centre of the main halo. This effect is apparent in the spin pa- 
rameter distribution which matches that of field haloes at larger 
radii but has a broader width than other published fits, extending 
to lower spin values. This difference between the spin properties of 
subhaloes and field haloes needs to be taken into account if precise 
measurement of the spin parameter distribution are to be made. 

The recovered spin parameter goes through a minimum for 
subhaloes near the edge of the host at about half the r,nax value. 
Here, if outer particles are stripped tidally as a substructure falls 
into a host halo, the result will be a decrease in the spin. This im- 
plies a radial dependant factor needs to be taken into account when 
compiling substructure catalogues, as the infalling haloes tend to 
have their outer particles removed. Once the outer layer has been 
lost the spin parameter generally increases to smaller radii as less 
and less mass is considered. 

The value of the spin parameter measured is dependent upon 
the choice of where to place the outer edge and precisely which 
material is included in the calculation. As we have shown here 
and elsewhere these choices are very halo finder dependent and so 
care should be taken when inter-comparing spin parameter mea- 
surements from different codes. 

In a future project we plan to look more closely at the differ- 
ence between field and substructure haloes, to compare more di- 
rectly the spin parameter found. 
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